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GENERALIZATION 

OF THE VAN DER WAALS EQUATION 
PACS 61.43.Gt, 64.70.F FOR ANISOTROPIC FLUIDS IN POROUS MEDIA 


The generalized van der Waals equation of state for anisotropic liquids in porous media con¬ 
sists of two terms. One of them is based on the equation of state for hard spherocylinders in 
random porous media obtained from the scaled particle theory. The second term is expressed 
in terms of the mean value of attractive intermolecular interaetions. The obtained equation 
is used for the investigation of the gas-liquid-nematic phase behavior of a molecular system 
depending on the anisotropy of molecule shapes, anisotropy of attractive intermolecular inter¬ 
actions, and porosity of a porous medium. It is shown that the anisotropic phase is formed 
by the anisotropy of attractive intermolecular interactions and by the anisotropy of molecular 
shapes. The anisotropy of molecular shapes shifts the phase diagram to lower densities and 
higher temperatures. The anisotropy of attractive interactions widens significantly the coexis¬ 
tence region between the isotropic and anisotropic phases and shifts it to the region of lower 
densities and higher temperatures. It is shown that, for sufficiently long spherocylinders, the 
Uquid-gas transition is localized completely within the nematic region. For all the considered 
cases, the decrease of the porosity shifts the phase diagram to the region of lower densities and 
lower temperatures. 

Keywords: fluids in random porous media, gas-liquid-nematic phase transitions, van der 
Waals equation, scaled particle theory, hard spherocylinders. 


1. Introduction 

It is a great pleasure for us to present our paper for 
publication in this special issue dedicated to the 70- 
th anniversary of Academician L.A. Bulavin, a known 
Ukrainian scientist in the physics of liquid state. His 
contribution to the development of the physics of liq¬ 
uid state is important. Some results obtained by him 
are partially summarized in books [ll-@| . At the same 
time, the progress of his former students inside and 
outside of Ukraine illustrates the importance of the 
scientific school created by L.A. Bulavin. 

The first understanding of the nature of the liq¬ 
uid state of matter is connected with the van der 
Waals equation formulated nearly 150 years ago 
This equation provided the possibility to describe 
the phase transition from the gaseous to the liquid 
state and to account for the presence of the criti¬ 
cal point beyond which the gaseous phase can not 
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be transformed into a liquid. It also provided the 
possibility to describe the coexistence between the 
liquid and gaseous phases, and to predict the ex¬ 
istence of metastable states, namely a supercooled 
gas and a superheated liquid. The background of the 
van der Waals equation is based on the idea of dif¬ 
ferent treatments of short-range repulsive and long- 
range attractive intermolecular interactions. The re¬ 
pulsive interactions fix the size and the shape of 
molecules and essentially determine the structural 
and entropic properties. The contribution of attrac¬ 
tive interactions is mainly energetic and can be de¬ 
scribed in the framework of the mean field approxi¬ 
mation. The first strong statistical mechanics treat¬ 
ment of the van der Waals equation was done about 
50 years ago for the hard sphere model with attrac¬ 
tive interactions in the form of the Kac potential [3,1^ 




( 1 ) 
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where r is the iiiterparticle distance. 

In the limit 7 —^ 0 , the equation of state can be 
presented in the form 0i| 

-na,,. 

P \p Ms 


( 2 ) 


where /3 = \/{kT), k is the Boltzmann constant, T is 
the temperature, P is the pressure of the fluid, p is the 
density, is the contribution of hard spheres 

(HS), T] = ^irpD^ is the fluid packing fraction, and 
D is the diameter of a hard sphere. The second term 
in Eq. ([5]) describes the contribution from attractive 
interactions through the constant a: 


a = — 


1 


LXJ 

J 7^C/(7r)r^dr. 


(3) 


D 


The equation of state in the form ([2|) is a gener¬ 
alization of the van der Waals equation. It coincides 
with it exactly in the one-dimensional case where the 
contribution from hard spheres is described by the 
Tonks equation 0- However, we can use a more cor¬ 
rect description of the hard sphere contribution such 
as, for example, the Carnahan-Starling equation |ll|. 
The equation of state in the form © can be used 
also for the description of non-spherical molecules and 
can be generalized for the description of isotropic and 
anisotropic fluids in porous media. This is the aim of 
this paper. 

To this end, we will use the Maddeii-Glandt model 
[T3 |. According to this model, a porous medium is 
presented as a quenched configuration of randomly 
distributed hard spheres forming a matrix, in the free 
space of which there are fluid molecules. A specific 
description of a fluid in such porous media is related 
to double quenched-annealed averages: the annealed 
average is taken over all the fluid configurations and 
an additional quenched average should be taken over 
all realizations of the matrix. 

The analytical results for a hard sphere fluid in 
hard sphere matrices 15,1^ obtained recently by ex¬ 
tending the scaled particle theory (SPT) pro- 

yide a strong basis for a generalization of the yan 
der Waals equation for simple fluids in porous nie- 
dia ll^. T he generalization of the results obtained 
in Il5l.ll(lll t o non-spherical molecules in porous me¬ 
dia im allowed us to generalize the yan der 
Waals equation [13 to anisotropic fluids in porous 
media. The inyestigations of the gas-liquid-nematic 


phase equilibria in the framework of the generalized 
yan der Waals equation show a rich yariety of phase 
behayiors that depends on the molecular shape, yalue 
of attractiye intermolecular interactions, and poros¬ 
ity of porous media. In this paper, we will continue 
this inyestigation. We will focus on the role of the 
anisotropy of attractiye interactions. 

The paper has the following structure. In Section [21 
we present the results for a hard spherocylinder fluid 
in porous media. In Section |3l we use these results to 
generalize the yan der Waals equation. In Section 21 
we study the influence of the anisotropy of intermolec¬ 
ular interactions on the phase behayior of a molecular 
fluid in a porous medium, by using the generalized 
yan der Waals equation. 


2. Application of the Scaled 
Particle Theory to the Description 
of Thermodynamic Properties 
of a Spherocylinder Fluid 
in a Random Porous Medium 


A hard spherocylinder fluid is widely used for the 
description of the influence of the molecular shape 
on the orientational ordering in anisotropic fluids 
0 ,( 23 . In this section, we apply the scaled parti¬ 
cle theory to the description of the thermodynamic 
properties of hard spherocylinders in random porous 
media created by hard spheres. The key point of the 
SPT theory is based on the deriyation of the chem¬ 
ical potential of an additional scaled particle of a 
yariable size inserted into a fluid. This excess chem¬ 
ical potential is equal to the work needed to cre¬ 
ate a cayity in a fluid, which is free from any other 
particles. The theory combines the exact considera¬ 
tion of an infinitely small particle with the thermo¬ 
dynamic consideration of a scaled particle of a suffi¬ 
ciently large size. The exact result for a point scaled 
particle in a hard sphere fluid confined in a random 
porous medium was obtained in 21|. Howeyer, this 
approach named SPTl contains a subtle inconsis¬ 
tency appearing, when the size of matrix particles 
is quite large compared to the size of fluid parti¬ 
cles. Later on, this inconsistency was eliminated in 
a new approach known as SPT2 0 | . In this section, 
we generalize this approach for hard spherocylinder 
fluids in random porous media. 

Following [l^[^.l^. we introduce an additional 
spherocylinder with the scaling diameter Ds and the 
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scaling length Lg into a spherocylinder fluid in a 
porous medium: 

Dg = XgDi^ Lg — o^gLi^ (4) 

where Di and Li are, respectively, the diameter and 
the length of the fluid spherocylinder. The excess 
chemical potential for a small scaled particle can be 

written in the form [13 


where P is the pressure of the fluid, and Vg is 
the volume of the scaled particle. The multiplier 
1 /po(q^s) As) appears due to an excluded volume occu¬ 
pied by matrix particles. The probability po(c^s> As) is 
directly related to two different types of porosity in¬ 
troduced by us in The first one corresponds 

to the geometrical porosity 


00 = Po (as = As = 0), 


( 9 ) 




lnpo(as, As) - In 


1 - 


m 


ViPo(as, As) 


characterizing the free volume for a fluid. For a hard 
sphere matrix. 


X (—-|- As)^ -l- —DiLi(l + As)^(l -l- as) -I- 


+ + As)as X 


X //(fli)/(n2) sin'di2(inidn2 


( 5 ) 


where rji = piVi is the fluid packing fraction, pi is the 
fluid density, Vi is the spherocylinder volume, and 

Po(as, As) = exp[-/3/i°(as, As)] (6) 

is the probability to And a cavity created by a scaled 
particle in the empty matrix. It is defined by the ex¬ 
cess chemical potential p,^(as, As) of the scaled parti¬ 
cle in the limit of an infinite dilution, fl = (d, p) de¬ 
notes the orientation of particles and is defined by the 
angles and (p; dfl = sin ■dd'ddp is the normalized 
angle element, 1912 is the angle between orientational 
vectors of two molecules, and /(fl) is the singlet ori¬ 
entational distribution function normalized in such a 
way that 

J f{n)dn = i. (7) 


<('0 = 1- Vo: (10) 

where ryo = is the packing fraction of the 

matrix, po is the density of matrix particles, and Dq 
is the diameter of matrix particles. 

The second type of porosity corresponds to the case 
As = as = 1 and leads to the thermodynamic porosity 

0 =Po(as = As = 1) = exp(-0/r?) (11) 

defined by the excess chemical potential of fluid parti¬ 
cles p)* in the limit of infinite dilution. It characterizes 
the adsorption of a fluid in the empty matrix. In the 
case under consideration. 


0 = (1 - ryo) exp 
9 


Vo 


I-Vo 


'T o('7l + 1) + 371 "T - 


vl ^ 2 

2971 ^ 


(1 - VO? 2 
(^(37r-l)V(l + po+p„^) 


( 12 ) 


where t = and 71 = 1 -I- Li/Hi. 

In accordance with the ansatz of the SPT theory 
[l3l - [l3 - w{Xg,ag) can be presented in the form 


We note that, here and below, we use conventional 
notations [IsMlSlj . where the index “1” is used to de¬ 
note a fluid component, the index “0” denotes matrix 
particles, while, for the scaled particles, the index “s” 
is used. 

For a large scaled particle, the excess chemical po¬ 
tential is given by the thermodynamic expression for 
the work needed to create a macroscopic cavity inside 
a fluid and can be presented in the form 


w{Xg,ag) = Woo + wioXg + 


+ WQiUg + WiiUgXg + 


W2oA^ 

2 


(13) 


The coefRcients of this expansion can be found from 
the continuity of and the corresponding deriva¬ 
tives dpf^/dXg, dpf^/dag, d^p.f^/{dag){dXg) and 
d'^pf^/dXl at As = Qfs = 0. After setting As = as = 1, 
we found the relation between the pressure P and the 
excess chemical potential of a fluid: 


PpT = w{as,Xg) + pPVg/poiXg, as), (8) /3 (pf - p°) = - In (1 - ryi/0o) -f x 
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+ (14) 


1 - 
where 

MHf)) = 6 + 


(l-??i/^o)^ Pi </> 


6 (71 - 1) t(/) 

371-1 


^ (4 3 ( 71 -) r(/) \ _ ^ / 671 

4>0 \ 371-1 j (j)o \ 371 -1 

p'dax _ 1 p'dw I n PoaPoX fP'ox 

(j)o 2 00 00 


13(r(/)) = 


671 


p'A/3 (271-1) 


v 37 i - 1 00 y y 371 - 1 

3 (71 - 1)^ T{f) p'o^ _ IppA 

371 - 1 00 2 00 /’ 


(15) 


(16) 




dpijr^ \dpijj,’ 

one derives the fluid compressibility in the form 


(18) 


^ \dpijj, (1- 771 / 0 ) 
+ (l + yl(r(/))) 


Vi/(t>o 


(1- 771 / 0 ) (I- 771 / 0 O) 


+ (^(r(/))+ 2 i?(r(/))) 


{vi/M 


(1 - 771 / 0 ) (1 - 77l/0o)' 


- 2 i?(r(/)) 


irii/(PoT 


(1 - 771 / 0 ) (1 - 77l/0o) 


(19) 


— 00 1 — 77l/0o 
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( 44 (r(/))+ 2 B(r(/)))- 


0—00 

^2 


m/'/o 

4> — 4>o Vi — Vi /'/'o 
In , 1 - 771 / 0 ^ + 2 g(r(/))x 


0 ??l/0O 

0 — 00 1 — 77l/0o 


4> — 4‘o 1 —^i/'/o 

'1 (??i/</o)^ 


2 (1 - 77i/0o)2 
1 - 771 / 0 ' 


■ In 


(0-0o)2 I- 771 / 0 O 

pp ^ l-?7i/</> 


( 20 ) 


Pi 

X 


m 


■ In 


1 - Pl/00 
1 - 771/0 


(1 + A(t(/))) X 


771 0 - 00 1 - 77l/0o 

(A(r(/))+ 2 B(r(/)))- 


r(/) = /(ni)/(fl2)sindi2dflidfl2- ( 17 ) 

/ _ Opo(as,Aa) / _ dpo{as,Xs) ri _ d'^PoiasAa) 

Pox — dXs ’ ^Oct — da^ ’ -^OaA — 9a«9As ’ 
^ Poias,Xs) corresponding derivatives 

at as = Xs = 0. 

Using the Gibbs-Duhem equation, which relates 
the pressure P of a fluid to its total chemical potential 
Pi = ln(Af Alii) + pf^, 


0-00 [1 — 77l/0o 


771 0 — 00 
■ 2 P(r(/)) 
1 


■In 


1 - 771/0 


1 - Pl/00 
1 


0-00 


m/Po 


20 - 00 


2(1-77i/0o)2 0-, 

1 - pl/</ 


■In 


l-pi/00 77 i(0-0o) 2 l-77l/0o_ 

where 

a(/) = J /(U)ln/(U)dU. 


( 21 ) 


( 22 ) 


As noted in 15 , expressions (I20h - (I21I1 have diver¬ 

gences at 771 = 0 and pi = 0o. Since 0 < 0o, the 
divergence at 771 = 0 occurs at lower densities and 
should be removed. Different corrections improving 
the SPT2 results were proposed in In this 

paper, we consider only the SPT2b approach, which 
can be derived if 0 is replaced by 0 o everywhere in 
HU) except the first term. In consequence, the chem¬ 
ical potential and the pressure of the fluid can be 
presented in the form 

/3(pr - = cr{f) - ln(l - 771 / 0 ) -h 

Pi/</o 


where Ai is the fluid thermal wave length, and 
A/j^ is the rotational partition function of a single 
molecule ( 2 ^. 

After the integration of relation m over pi, one 
obtains the expressions for the chemical potential and 
for the pressure in the SPT2 approximation [iM3: 

PipT - Pi) = - ln(l - Pi/0 ) + (^(t(/)) + 1 ) X 

1- 771/0 


+ (1 + A(r(/))) 


1 - Pl/00 


-b 


+ l(A(.(/)) + 2B(r,/)))U5!^ 


(23) 


■/3PY 


SPT2b 


Pi 


00 


m 


PiJ Pi V 9/ Pi V 90 

1 , ^(t(/)) Pi/4'0 

l-pi/00 2 (1 - 77 i/0o)2 
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2B{T{f)) fa/0o)^ 

3 (1 - 


(24) 


which reproduces quite well the results of computer 
simulations [l5l - [l8l . f2H | . 

From the thermodynamic relation 

0F 

— = Puipi - I3P, (25) 

one can obtain the expression for the free energy: 


^y-ij^SPT 2 b ^ ^ Pi(ln(A?pi) - 1) + 


+ /3MiPi - Pi ln(l - ?7i/</>) + 

+ — ln(l - rii/4>) - ln(l - rji/M ■ 
Vi m 


^('r(/)) Vi/(po 
^ 2 l-77i/(/)o 

3 \l-r]i/ 4 'o)' 


(26) 


By minimizing the free energy with respect to the 
variation of the distribution function /(S 2 ), we obtain 
the integral equation 


ln/(Oi) + l + C 


f{^2) sindi2(in2 = 0, 


(27) 


In the presence of a porous medium within the On- 
sager model, we have 

c,/0o = 3.290, cn/^o = 4.191. (31) 

For finite Li and Di, 

a = 3.290, Cn = 4.191, (32) 

where Ci and Cn are defined by (l2^ . These values 
of Ci and Cn define the isotropic-nematic phase dia¬ 
gram depending on the ratio Li/Di and the matrix 
parameters for a hard spherocylinder fluid in a ma¬ 
trix. As was shown in [I 3 , the obtained theoretical 
results are in agreement with the data of computer 
simulations (^ . 

3. Generalization of the van der Waals 
Equation for Anisotropic Fluids in Random 
Porous Media 

We will use the results for a hard spherocylinder fluid 
presented in the previous section as a reference system 
for the generalization of the van der Waals equation 
for anisotropic fluids in random porous media. 

Such generalization includes the non-spherical 
shape of molecules, anisotropy of the intermolecular 
interaction, and presence of a porous medium. As a 
result, we will have a more general form of Eq. m 

[13: 


where 

^ ^ ml( 1^0 \ 3(71 - 1 )^ A _ Pox \ , 

1 - m/(l>o [ 371-1 V 2 (j)oJ 
m/ 4'0 ill - 1)^ f 671 _ p'g^\ 
(1 - m/ 4 o) 371 -1 V371 -1 4 o)_ 


(28) 


This equation can be solved numerically, by us¬ 
ing an iteration procedure according to the algorithm 
proposed in [^. We note that, in the Onsager limit 
Li — 00 , Di ^ 0 [l9| . 

C —^ c = —ttT^ZIiPi, (29) 


where c is finite. 

From the bifurcation analysis, it is found that 
Eq. ((Tfl) has two characteristic points [53: 


Ci = 3.290, Cn = 4.191, 


(30) 


where Ci corresponds to the highest density of a sta¬ 
ble isotropic fluid, and c„ is related to the minimum 
density of a stable nematic fluid. 


Pi 


—) - 12?7ia/3, 

. Pi /HSC 


(33) 


where (—) is the contribution from hard sphero- 

cylinders (HSC) in porous media, which is described 
by expression (IMl) . The contribution from attractive 
interactions is described by a constant a, which can 
be presented in the form 


1 

Wi 


f f{n,)f{n2)u^^\n2n,n2)dm2dn,dn2, 

(34) 


where the factor ^ excludes the volume occupied by 
matrix particles, Vi is the volume of a molecule, m = 
piVi, and C/®'**‘'(ri 2 Hifl 2 ) is the attractive part of the 
intermolecular interaction. 

Similarly to [i3, we present the potential 
{r 12^10,2) in the form of a modified Lennard- 
Jones potential 


C/"“(ri 2 Hi 02 ) 


C/lj 


fj(HiH 2 iij.) 


P 12 


X 
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X [1 +X-P 2 (cOS 1 ? 12 )] , 


(35) 


C/lj = 4eo 


ri2 


r(r 2 ir 22 rir) 

ri 2 


(7 ( ] i ^2 

r \2 


12 


, ri2 > tT(r2ir22f^r), 


^ ^ ^r> n r> 1 

(7lj I-) = 0, ri2 < Cr(r2lS22iir), 


?'12 


(36) 


where P 2 (cosi?i 2 ) = 5 ( 3008^1912 — 1) is the sec¬ 
ond Legendre polynomial, t9i2 is the angle be¬ 
tween the principal axes of two interacting molecules, 
and CT(nin2f9r) is the contact distance between 
molecules. It depends on the orientations of two 
molecules, fli and 172, as well as on the orientation 
of the distance vector fi 2 between their centers. In 
the case where the repulsive part of the interaction 
is spherical {Di = u), this potential reduces to the 
Maier-Saupe potential [ 2 ^. We note that potential 
(I5S)) is the sum of two Lennard-Jones potentials. The 
first one is related to the isotropic attraction, and the 
second one to the anisotropic attraction. The ratio 
of the well depths of these two potentials y = £ 2/60 
specifies the degree of anisotropy in the attraction of 
the total potential. 

Following the traditional scheme [s^, taking into 
account that df = r'^drdflr, and using a dimension¬ 
less intermolecular distance r* = r/a{nift 2 ^r), one 
obtains 


1 


df7idf72/(I7i)/(I72) [1 + xP2(cosi9i2)] x 


X F“^(I7iI72)3 J r*^dr*pUu{r*), (37) 

0 


where 

I4xc(f^lf92) = - 


dI7r[o'(I7iI72^^r)]^ 


(38) 


is the excluded volume formed by two spherocylinders 
with orientations I7i and 172- 

The following expressions for the chemical potential 
and the free energy correspond to Eq. (|55)) : 


- d-i) = cr{f) - ln(l - r]i/(j)) + 




— = Pict(/) + pi(ln(A?pi) - 1) -f - 

- Pi ln(l - r]i/(j)) + — ln(l - pi/fj)) - 

m 


Vi 


ln(l - rii/(j)o) + pi 


2 l-?7i/(^o 


-pi- 




f vi/'/>o Y 
Vi - m/w 


12piapi/3. 


(40) 


From the last expression, we have the following inte¬ 
gral equation for the singlet distribution function: 


In/(I7i)-|-1-|-C J /(I72) sin i9i2dI72-f 

-I- J /(I72)P2(cOSl9i2)df72 X 

• f a(I7iI72l7,)7 ^ ^ 


PpiX 


(41) 


where C is given by expression (l2^ . 

The obtained singlet distribution function /(f7) is 
used in (IM)) for the calculation of the parameter a. 


4. Influence of Anisotropic 
Attractive Intermolecular Interactions 
on the Phase Behavior of Anisotropic 
Fluids in Porous Media 

Now, we apply the developed theory to the descrip¬ 
tion of the gas-liquid-nematic phase behavior of the 
considered molecular fluids in porous media created 
by a random configuration of hard spheres. Given the 
chemical potential and the pressure as functions of 
the density pi at different temperatures, one can cal¬ 
culate the coexistence curves from the conditions of 
thermodynamic equilibrium: 

di{p\.T) = pi{pl,T), 

P{p\,T)=P{plT), 

where p\ and p\ are the fluid densities of two different 
phases 1 and 2. The numerical solution of these equa¬ 
tions is realized with the use of the Newton-Raphson 
algorithm. 
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In contrast to (l^ . where our investigations were 
concentrated on the influence of the molecular shape 
on the phase behavior of molecular fluids with x = 0, 
we will focus more on the influence of the anisotropy 
of attractive intermolecular interactions. The corre¬ 
sponding results of our investigation of the influence 
of parameter x different values of porosity (jjQ and 
different values of parameter Li/Di are presented in 
Figs. [THU The results are presented in the form of the 
phase diagram in the coordinates “dimensionless den¬ 
sity T]i - dimensionless temperature T* = kT / eg”. In 
order to compare the influence of the parameter x, we 
also present the results for y = 0. We note that, in 
accordance with dSZl), the contribution of anisotropic 
attractive interactions is proportional to 

y/(ni)/(fl2)^2(cosdi2)dfIidfl2 = Sf (43) 

where 

S2 = J P2(cosd)/(fI)dfI (44) 

is the order parameter. 

Since S '2 = 0 in the isotropic phase, the influence of 
anisotropic attractions in the isotropic phase is negli¬ 
gible in the van der Waals approximation. In order to 
check the law of corresponding states in each figure, 
we present the phase diagrams also in the reduced 
variables rji/rjic — T* /T*, where rjic and T* are the 
corresponding values of critical density and critical 
temperature for the gas-liquid phase transition. 

We start from the hard sphere model with attrac¬ 
tive interactions in the form ((36l) . In this case, L = 0, 
and D = a is the diameter of hard spheres. The cor¬ 
responding phase diagram is presented in Fig. [TJ In 
the region of small densities, we have the gas phase G, 
which changes, as the density increases, to the liquid 
phase L and, at higher densities, to the anisotropic 
nematic phase N. The nematic phase appears due to 
the anisotropy of attractive interactions, and the bi¬ 
furcation line of the nematic phase is proportional 
to , according to m- This means that the 
phase transition appears at higher densities and vice 
versa, as the temperature increases. As the tempera¬ 
ture decreases, the phase transition appears at lower 
densities. At sufficiently low temperatures, the region 
of the gas-liquid transition converges to the tricriti- 
cal gas-liquid-nematic {G — L — N) point. Below the 



a 



b 

Fig. 1. Temperature-density phase diagram calculated from 
the generalized van der Waals equation for hard spheres with 
isotropic and anisotropic attractive interactions (x = 0.3) in 
random porous media with the porosity = 0.8 and r = 
= 1/3. For the purpose of comparison, the phase diagram in 
the absence of a porous medium 0o = 1.0 is also presented. The 
dashed lines correspond to the case x = 0. On the top (case 
a), the phase diagram in the coordinates T* = kT/eo and 
rji = piVi is presented. At the bottom (case b), the phase 
diagram is presented in the coordinates T* jT* and r\\lr\\^c 

tricritical point, only the gas-nematic (GN) coexis¬ 
tence is seen. At temperatures higher than the tricrit¬ 
ical one, the anisotropic attractive interaction does 
not change the gas-liquid coexistence line in the van 
der Waals approximation. But, for lower tempera¬ 
tures, the anisotropic attraction slightly widens the 
gas-liquid coexistence and leads to a strong widen- 
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a 



ll/Tll.c 

b 

Fig. 2. Temperature-density phase diagram calculated from 
the generalized van der Waals equation for hard spherocylin- 
ders with the anisotropy of sizes Li/Di = 3 and with isotropic 
and anisotropic attractive interactions (x = 0.3) in random 
porous media with the porosity (f>o = 0.8 and r = 1/6. For the 
purpose of comparison, the phase diagram in the absence of a 
porous medium (f>o = 1.0 is also presented. The dashed lines 
correspond to the case x = 0. On the top (case a), the phase 
diagram in the coordinates T* = kTjeo and rji = piVi is pre¬ 
sented. At the bottom (case b), the phase diagram is presented 
in the coordinates T* jT* and 


ing for larger densities. The presence of a porous 
medium shifts the phase diagram to the region of 
lower densities and lower temperature, as the poros¬ 
ity 00 decreases. In contrast to the usual gas-liquid 
phase transition ^ , the law of corresponding states 


in the considered case is more or less valid in the 
region of small densities. But, at higher densities, 
the coexistence curve becomes wider, as the poros¬ 
ity decreases. 

The influence of the parameter L\/D\ responsible 
for the anisotropy of a molecular shape is presented 
in Figs.[2Hll Comparing Fig.[l]and Fig. [21 we can see 
that the non-sphericity of a molecular shape leads 
to a shift of the phase diagram to lower densities and 
higher temperatures. The increase of the ratio Li/Di 
leads to an increase of the critical temperature and 
a decrease of the critical density for the gas-liquid 
transition. The tricritical temperature increases also, 
and the gas-liquid region is essentially narrower than 
in the case of spherical molecules. Similarly as for 
Li/Di = 0 (Fig. [T|), the anisotropic phase appears 
in the case Li/Di = 3 (Fig. |2]) due to the anisotropy 
of attractive intermolecular interactions. However, 
the anisotropy of molecule shapes highly modifies 
the region of the coexistence of isotropic and ne¬ 
matic phases. The liquid-nematic coexistence region 
becomes wider and is not very sensitive to the tem¬ 
perature. Similarly as for Li/Di = 0 below the tri¬ 
critical temperature, the anisotropic attractive in¬ 
teraction does not change the gas-liquid coexistence 
line. For lower temperatures, the coexistence region 
slightly widens at lower densities and widens signifi¬ 
cantly at higher densities. We note that, at x = 0, the 
nematic phase does not appear. I n ag reement with 
the data of computer simulations (32|, the nematic 
phase in the fluid of hard spherocylinders appears 
at Li/Di > 3.7. For L\/D\ = 3 similarly as for 
Li/Di = 0, the decrease of the porosity 0o shifts 
the phase diagram to lower densities and lower tem¬ 
peratures. In contrast to the case Li/Di = 0, the 
corresponding law in the case Li/Di = 3 is more or 
less satisfied for all densities, including the region of 
the isotropic-nematic transition. 

The phase diagram for the case Li/Di = 5 is pre¬ 
sented in Fig. m Such asymmetry of the shape of 
molecules is sufficient for the nematic ordering. As 
we noted in [17|, the anisotropic attractive interac¬ 
tion expands the region of orientational ordering. The 
liquid-nematic region widens significantly if x iH” 
creases. As a consequence, the gas-liquid transition 
disappears for a sufficiently large anisotropy, and only 
the isotropic-nematic transition takes place. Such sit¬ 
uation for Li/Di = 5 is observed at x = 0.3. Due 
to this fact, we put x = O.f in Figs. |3|and|Tl As we 
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Fig. 3. Temperature-density phase diagram calculated from 
the generalized van der Waals equation for hard spherocylin- 
ders with the anisotropy of sizes LijDi = 5 and with isotropic 
and anisotropic attractive interactions (x — 0-3) in random 
porous media with the porosity 0o = 0.8 and r = 0.1. For 
the purpose of comparison, the phase diagram in the absence 
of a porous medium, (j!)o = 1.0, is also presented. The dashed 
lines correspond to the case y = 0. On top (case a), the phase 
diagram in the coordinates T* = kT/eo and 771 = piVi is pre¬ 
sented. At the bottom (case b), the phase diagram is presented 
in the coordinates T* jT* and 

can see from Fig. |3l similarly to Figs. [1] and [2l the 
gas-liquid coexistence line below the tricritical tem¬ 
perature does not depend on the anisotropic attrac¬ 
tive interaction. The anisotropic attractive interac¬ 
tion leads to an increase of the tricritical temperature 

ISSN 2071-0186. Ukr. J. Phys. 2015. Vol. 60, No. 8 


b 

Fig. 4' Temperature-density phase diagram calculated from 
the generalized van der Waals equation for hard spheres 
LijDi = 80 with isotropic and anisotropic attractive inter¬ 
actions (x = 0.3) in random porous media with the porosity 
(pQ = 0.8 and r = 0.01. For the purpose of comparison, the 
phase diagram in the absence of a porous medium, (po = 1 . 0 , 
is also presented. The dashed lines correspond to the case 
X = 0. On top (case a), the phase diagram in the coordinates 
T* = kT/eo and 771 = piVi is presented. At the bottom (case 
b), the phase diagram is presented in the coordinates T* jT* 
and riilni^c 

and a decrease of the tricritical density. As a result, 
the gas-liquid region is narrower, while the liquid- 
nematic region covers a wide range of densities. The 
presence of a porous medium shifts the phase diagram 
to lower densities and lower temperatures. The cor- 
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responding law is satisfied rather well. In the region 
of the anisotropic transition, the corresponding law is 
satisfied separately for y = 0 and for y = 0.1. 

Finally, we consider Fig. |31 where the phase dia¬ 
gram for Li/Di = 80 is presented. As we noted in 
|l7| . the transition into the nematic phase shifts in 
this case at y = 0 to the region of small densities. As a 
consequence, the gas-liquid transition appears in the 
nematic region. In accordance with the classification 
in , the gas phase in the nematic region is marked 
as Nj, and the liquid phase in the nematic region is 
marked as Nn. In contrast to the previous figures, 
in the case considered in Fig. 01 the influence of an 
anisotropic attractive interaction is of importance for 
the entire phase diagram. The anisotropic attractive 
interaction does not change significantly the critical 
density of the gas-liquid transition, but it changes 
strongly the value of critical temperature. In contrast 
to the case Li/Di = 5, the anisotropic attractive in¬ 
teraction in the considered case Li/Di = 80 leads to 
a strong increase of the tricritical temperature and 
to the expansion of the gas-liquid coexistence. The 
region of coexistence of the isotropic and nematic 
phases is also expanded and shifts to lower densi¬ 
ties. In the presence of a porous medium, the phase 
diagram is shifted to lower densities and lower tem¬ 
peratures, as the porosity (jjQ decreases. The corre¬ 
sponding law is satisfied quite well for the isotropic 
phase independently of the value of x- However, the 
values of y and porosity (po have a strong influence 
on the phase diagram in the anisotropic region. At 
the critical point, the corresponding law is satisfied, 
but, for higher densities, the behavior of the phase 
diagram depends significantly on the value of poros¬ 
ity po- 

5. Conclusion 

In this paper, we have presented the generalized van 
der Waals equation for anisotropic molecular fluids 
in porous media. This generalization is based on an¬ 
alytical expressions for the equation of state and the 
chemical potential of a hard spherocylinder fluid in 
a random porous medium obtained in the framework 
of the scaled particle theory. The second term of the 
generalized van der Waals equation is the mean value 
of attractive intermolecular interactions. By minimiz¬ 
ing the free energy of the fluid, we have obtained a 
nonlinear integral equation for the singlet distribution 


function, which describes the orientational ordering in 
the system. This ordering is connected with the non¬ 
sphericity of molecular shapes and the anisotropy of 
intermolecular attractive interactions. 

The investigations based on the generalized van der 
Waals equation demonstrate a wide variety of gas- 
liquid-nematic phase behaviors in molecular systems 
depending on the anisotropy of a shape of molecules, 
anisotropy of attractive interparticle interactions, and 
porosity of a porous medium. Here, we have focused 
our attention on the influence of the anisotropy of at¬ 
tractive interactions, which is the main reason for the 
orientational ordering in quasispherical molecules. It 
is shown that the anisotropy of molecular shapes 
leads to a shift of the critical point of the gas-liquid 
transition to lower densities and higher tempera¬ 
tures and to an increase in the tricritical temper¬ 
ature. The anisotropy of molecular shapes also ex¬ 
pands significantly the region of coexistence of the 
isotropic and nematic phases. It is shown that, for a 
larger anisotropy, the nematic phase can appear due 
to the anisotropy of molecular shapes. In this case, 
the anisotropy of attractive intermolecular interac¬ 
tions expands significantly the region of coexistence 
between the isotropic and nematic phases and shifts 
it to the region of lower densities and higher tem¬ 
peratures. Finally, at a sufficiently large anisotropy 
of molecular shapes, the transition into the nematic 
phase takes place at very low densities. As a result, 
the gas-liquid transition takes place in the nematic 
region. In all the considered cases, the decrease of 
the porosity of a porous medium shifts the phase di¬ 
agram to the region of lower temperatures and lower 
densities. 

The authors express their gratitude to Taras Pat- 
sahan and Ivan Kravtsiv for useful discussions and 
important suggestions during the preparation of this 
paper. 
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ySAFAJIBHEHHH PIBHHHHH 
BAH-/];EP-BAAJIbCA HA AHISOTPOHHI 
Pi;i,HHH B HOPHCTHX CEPEnORHTHAX 

P e 3 lo M e 

Hpe^CTaBjiene ysarajibHene piBHHHHH Ban-flep-Baajibca na 
aniaoTponHi pi^HHH b nopncxHx cepeflOBHU],ax CKjia;i;aeTbCH 3 
flBOx flO^aHKiB. HepmHH 3 hhx OaayeTbcn na piBHanni cxany 
TBep;i;Hx c4)epoLi,HjiiH^piB y BHna;i;KOBOMy nopncTOMy cepe;i;o- 
BUTTp, OTpHMaHHH B paMKax MBTO^y yaarajibHeHHH Macmxa- 
6Hoi" uacTHHKH. ^pyrHH ;i;o;i;aHOK BnpaiacaeTbCH uepea cepe^ne 
aHaueHHn noTenpiajiy npHTnrajibHoi’ MiiacMOjieKyjinpHoi" B3ae- 
Mo;i;h. Ha ocnoBi OTpHManoro piBHHHHH npoBe;],eHO ;i;ocjii;i;:>Ke- 
HHH 4)a30Boi" noBe^iHKH ra3—pi;i;HHa—HeMaTHK MOjieKyjinpnHx 
CHCxeM y 3ajie:acHOCTi Bi^ aHiaoTponh 4)OpMH MOjiexyji, aniao- 
Tponh npHTnrajibHoi" B3aeMo;i;h xa nopncxocxi nopncxoro cepe- 
^OBHU];a. HoKaaano, u],o aniaoxponna 4)a3a 4>opMyexbcn: hk 3a 
paxynoK aniaoxponHoi’ npHxnrajibHoi" B3aeMo;i;h, xax i 3a paxy- 
HOK aHiaoxponh 4)opMH MOjieKyji. Aniaoxponin 4)opMH MOjie- 
Kyji npHBo;i;HXb flo 3cyBy 4)a30Boi" ^iarpaMH b o6jiacxb Menmnx 
rycxHH xa bhttt,hx xeiunepaxyp, a aHiaoxponin npnxnrajibHoi" 
B3aeMo;i;ii anauHO poamnpioe o6jiacxb cniBicnyBaHnn iaoxpo- 
nnoi" xa neMaxHUHoi’ (^aa i xaKO»: acoBye ii b o6jiacxb hh:>khhx 
rycxHH i bhu];hx xeMnepaxyp. HoKaaaHO, mo npn ;i;ocxaxHbO 
BejiHKin acHMexph 4)OpMH MOjiexyji 4)a30BHH nepexm pmHHa— 
raa anaxo^HXbcn noBnicxio b o6jiacxi HeMaxHunoi" (^aan. Y Bcix 
BHna;i;Kax, mo poarjm^aioxbcn:, noHH:aceHHH nopncxocxi nopn- 
cxoro cepe;i;oBHm3' scoBye 4)a30By marpamy b o6jiacxb Hnncnnx 
rycxHH i xeMnepaxyp. 
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